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Abstract 

The program of systematic large-scale self-consistent nuclear mass calculations that 
is based on the nuclear density functional theory represents a rich scientific agenda 
that is closely aligned with the main research directions in modern nuclear structure 
and astrophysics, especially the radioactive nuclear beam physics. The quest for the 
microscopic understanding of the phenomenon of nuclear binding represents, in fact, 
a number of fundamental and crucial questions of the quantum many-body problem, 
including the proper treatment of correlations and dynamics in the presence of 
symmetry breaking. Recent advances and open problems in the field of nuclear 
mass calculations are presented and discussed. 
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1 Introduction 



The study of nuclei far from stability is an increasingly important part of 
a nuclear physics portfolio [1,2,3]. As radioactive beams gradually expand 
the borders of the nuclear landscape, theoretical modeling of the nucleus is 
changing in significant ways. The crucial question for the field [2], namely 
"What binds protons and neutrons into stable nuclei and rare isotopes?," 
nicely underlines this point: indeed, the data on rare isotopes with the large 
neutron-to-proton imbalance indicate that there are many gaps in our present 
understanding. 

Short-lived exotic nuclei offer unique tests of those aspects of the nuclear 
theory that depend on neutron excess [4,5]. The major challenge is to predict 
or describe in detail exotic new properties of nuclei far from the stability 
valley, and to explain the origins of these properties. New ideas and progress 
in computer technology have allowed nuclear theorists to understand bits and 
pieces of nuclear structure quantitatively. 

The new experimental developments inevitably require safe and reliable theo- 
retical predictions of nuclear properties throughout the whole nuclear chart in 
two main directions: (i) along the isospin axis, i.e., going outwards from the 
beta stability line to the neutron and proton drip lines, and (ii) towards the 
uncharted territory of super-heavy elements at the limit of mass and charge. 
The tool of choice is the nuclear density functional theory (DFT) based on 
the self-consistent Hartree-Fock-Bogoliubov (HFB) method. The key compo- 
nent is the universal energy density functional, which will be able to describe 
properties of finite nuclei as well as extended asymmetric nucleonic matter. 
The development of such a universal functional, including dynamical effects 
and symmetry restoration, is one of the main goals of the field. 

By employing various criteria (agreement with measured masses, radii, low- 
lying excited states, giant vibrations, rotational properties, and other global 
nuclear characteristics), one aims at adjusting the coupling constants of the 
functional. By finding correlations between parameters, one hopes to reduce 
their number and to understand physical reasons why different parametriza- 
tions yield similar results. One may also want to expand the parametrizations 
to cover aspects dictated by physics arguments and/or motivations coming 
from the effective field theory and QCD. The main challenges in this quest 
have been nicely summarized through five questions [6]: 

• What is the form of the nuclear energy density functional? 

• What are the constraints on the nuclear energy density functional? 

• What is the form of the pairing functional? 

• How to account for quantum correlations and symmetry-breaking effects? 
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• How to optimize computational techniques and error analysis? 

The aim of this paper is to briefly review the present state of the large-scale 
microscopic nuclear mass calculations and to discuss improvements needed. 
Section 2 introduces the DFT and Skyrme-HFB method. Some details con- 
cerning global mass calculations are given in Sec. 3. The long-term program 
is outlined in Sec. 4. Finally the summary is given in Sec. 5. 



2 Nuclear Energy Density Functional 

A theoretical framework aiming at the microscopic description of nuclear 
masses and capable of extrapolating into an unknown territory must fulfill 
several strict requirements. First, it must be general enough to be confidently 
applied to a region of the nuclear landscape whose properties are largely un- 
known. Second, it should be capable of handling symmetry-breaking effects 
resulting in a large variety of intrinsic nuclear deformations. Thirdly, it should 
describe finite nuclei and the bulk nuclear matter. Finally, in addition to ob- 
servables, the method should provide associated error bars. 

These requirements are met by the DFT in the formulation of Kohn and Sham 
[7]. The main ingredient of the non-relativistic nuclear DFT [8] (for relativistic 
nuclear DFT, see Ref. [9]) is the energy density functional that depends on 
densities and currents representing distributions of nucleonic matter, spins, 
momentum, and kinetic energy, as well as their derivatives (gradient terms). 
Standard Skyrme functionals employed in self-consistent mean-field calcula- 
tions are parametrized by means of about ten coupling constants that are 
adjusted to basic properties of nuclear matter (e.g., saturation density, bind- 
ing energy per nucleon) and to selected data on magic nuclei. The functionals 
are augmented by the pairing term which describes nuclear superfluidity [10]. 
When not corrected by additional phenomenological terms, standard function- 
als reproduce total binding energies with an rms error of the order of 2 to 4 
MeV [11,12,13]. However, they have been successfully tested over the whole 
nuclear chart to a broad range of phenomena, and usually perform quite well 
when applied to energy differences, radii, and nuclear moments and deforma- 
tions [8]. 

Historically, the first nuclear energy density functionals appeared in the con- 
text of Hartree-Fock (HF) or HFB methods and zero-range interactions such 
as the Skyrme force. However, it was realized afterwards that - in the spirit of 
the DFT - an effective interaction could be secondary to the functional, i.e., 
it is the density functional that defines the force. This is the strategy that we 
are going to follow. 
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2.1 The Densities 



The main ingredients of the nuclear DFT are the local nucleonic densities. 
Following the standard definitions [14,10], one considers local particle-hole (p- 
h) densities: particle p(r), kinetic r(r), spin Sk(r), spin-kinetic Tk{r), current 
jk(r), tensor-kinetic Fk(r), spin-current Jfcz(r), as well as the correspond- 
ing local particle-particle (p-p; or pairing) densities: p(r), f(r), Sk(v), Tfc(r), 
3k( r ), F k(r), and J fc ,(r). 

The local p-h and p-p densities are defined by the spin-dependent one-body 
density matrices: 

1 i 

p(r<r,rV) = -p(r,r')5 aa/ + - ^2(a\ai\a')pi(r, r'), 



1 1 

p(r<r,rV) = -p(r,r')5 aa/ + - ^{a\ai\a')pi{r, r'). 



For instance, 

p(r) = p(r,r), r(r) = V r V r 'p(r, r')\ r , = r , p(r) = p(r,r), 

(2) 

JuW = ^(V,-V9 Pj (r,r')| r=r . 

Since the nuclear DFT deals with two kinds of nucleons, the isospin degree of 
freedom has to be introduced and the isoscalar and isovector densities have to 
be considered [10]. 

2.2 The Energy Density Functional 

The energy density functional has the form 

E[p,p] = J d 3 rH(r), (3) 



where energy density TL{r) is usually written as a sum of the p-h energy density 
H(r) and the p-p energy density H(r). According to the DFT, there exists 
a nuclear universal energy functional that yields the exact binding energy of 
the nuclear system. This functional includes, in principle, all many-nucleon 
correlations. 



4 



The actual form of the nuclear energy functional is unknown. The strategy 
adopted by many practitioners is to build a functional around that generated 
by the Skyrme interaction. The most general form of the energy density func- 
tional that is quadratic in local densities and preserves the basic symmetries 
of the strong interaction, such as time-reversal symmetry, charge invariance, 
and proton-neutron symmetry, has been proposed in Ref. [10]. In practical 
applications, however, simplified forms of this functional have been used. For 
instance, one particular representation of the energy functional for the ground 
states of even-even nuclei can be written as:: 



+ \h [(l + \ Xl ) p (t - I Ap)] - (± + an) E Pq (r q - \A Pq ) 
+ \t 2 [(l + \x 2 ) p(r+ \Ap) + (| + x 2 ) E p q (r q + \A Pl 



+ ±t 3 ff \(l + \ XZ ) p 2 - (x 3 + i) E Pf\ - I (tlX! + t 2 X 2 ) E J 

y i j 

+ l(ti- 1 2 ) E JLj - \ Wo E % fc [pVfcJij + E P^fcJ,^-] 



(4) 



and 




(5) 



where q labels the neutron (q = n) or proton (q = p) densities and the 
quantities which do not carry index q are the isoscalar densities (sums of 
proton and neutron densities; e.g., p = p{r) = p n (r) + p p {r)). 

The p-p energy functional (5) corresponds to a density-dependent delta in- 
teraction. Usually, 7 = 1, p =0.16 fm -3 , and 14=0, 1, or 1/2 for volume-, 
surface-, or mixed-type pairing. In Eq. (4), H c (r) stands for the Coulomb 
energy density with the exchange term treated in the Slater approximation. 



As seen from Eqs. (4) and (5), typical Skyrme density functionals include 
about 14 unknown parameters. Some of them are usually adjusted to repro- 
duce the basic properties of the infinite nuclear matter while the remaining 
coupling constants are fitted to known nuclear masses, radii, and other mea- 
sured properties. 
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2.3 Variational Equations 



By varying the energy functional (3) with respect to the density matrices p 
and p one arrives at the HFB equations: 



= E \, (6) 





where U = U(E,ra), V = V(E,ra) are the HFB wave functions, and h and 
h are the local particle and pairing mean-field Hamiltonians. 

The HFB equations (6), also called the Bogoliubov de Gennes equations by 
condensed matter physicists, are the generalized Kohn Sham equations of the 
DFT. It is worth noting that - in its original formulation [15] - the DFT formal- 
ism implicitly includes the full correlation functional. In most nuclear applica- 
tions, however, the correlation corrections are added afterwards. Those correc- 
tions usually include the following terms: the center-of-mass correction, rota- 
tional correction associated with the spontaneous breaking of rotational sym- 
metry, vibrational correction (quantum zero-point vibrational fluctuations), 
particle-number correction due to the broken gauge invariance, as well as other 
terms. 

The spectrum of quasi-particle energies E is continuous for |i?|>— A and dis- 
crete for |i?|<— A. However, when solving the HFB equations on a coordinate- 
space lattice of points or by expanding quasi-particle wave functions in a finite 
basis, the quasi-particle spectrum is discretized and one can use the notation 
V k (ra) = V(E k ,ra) and U k (ra) = U(E k ,ra). Since for E k >0 and A<0 the 
lower components V k (ra) are localized functions of r, the density matrices, 



p(ra,r'a') = J2V k (ra)V k *(r'a'), p(ra,r'a f ) = - ]T V k (ra)U* k (r'a'), (7) 

k k 

are always localized. The norms N k of the lower components define the total 
number of particles 

N k = Y,Jd 3 r\V k (ra)\ 2 , N = £ N k = J d 3 r p(r). (8) 



For spherical nuclei, the self-consistent HFB equations are best solved in the 
coordinate space where they form a set of ID radial differential equations 
[16,17]. In the case of deformed nuclei, however, the solution of deformed HFB 
equations in coordinate space is a difficult and time-consuming task. For axial 
nuclei, the corresponding 2D differential equations can be solved by using the 
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basis-spline methods (see, e.g., Ref. [18]). For triaxial nuclei, 3D solutions in 
a restricted space are possible by using the so-called two- basis method [19]. 



3 Large-Scale Microscopic Nuclear Mass Calculations 

The large-scale microscopic nuclear mass calculations, such as those of Refs. 
[13,20,21,22], typically require that the variational equations are repeatedly 
solved for thousands of nuclei. For example, when adjusting the parameters 
of the energy density functional to measured masses, one has to calculate 
ground-state configurations of around two-thousand nuclei many times dur- 
ing the fitting process. Actually, the situation is even more complicated, as 




L=— i i i i i i i i i i i i i i i i i i l 

20 40 60 80 100 120 140 160 180 

Neutron Number 



Fig. 1. Quadrupole deformations (3 for all even-even particle-bound nuclei calculated 
with the SkP energy density functional [16] in the p-h channel and the volume delta 
pairing using the deformed HFB+THO code with 20 THO shells. 

several independent calculations have to be carried out for a given nucleus 
to find the ground-state energy of the system among several coexisting local 
minima. Furthermore, if odd-A and odd-odd nuclei are considered during the 
fitting process, many one-quasiparticle and two-quasiparticle states have to 
be considered to find the actual ground state. Finally, when the functional 
has been established, properties of around ten-thousand particle-bound nuclei 
throughout the nuclear chart can be computed. All in all, fitting a functional 
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and preparing a mass table is a challenging computational problem that re- 
quires highly optimized numerical codes and excellent utilization of modern 
multiprocessor computer resources. 




20 40 60 80 100 120 140 160 180 
Neutron Number 

Fig. 2. Similar to Fig. 1 except for two-neutron separation energies. 

Our group has been laying out theoretical foundations and constructing com- 
putational tools to tackle this ambitious task. We utilize a fast HFBRAD code 
for spherical HFB calculations [23], which takes no more than 10 CPU minutes 
per nucleus on an Intel Xeon 2.8 GHz processor, as well as the HFBTHO code 
for axially deformed HFB calculations [24,25] with acceptable processor speed 
- averaging to about 1 CPU hour per nucleus. 

The large-scale mass calculations based on the HFB+THO code, extended 
with a minimal MPI communication in order to run in a parallel regime across 
the nodes of the multiprocessor computer, are illustrated in Figs. 1 and 2, 
which display, respectively, calculated charts of nuclear deformations and two- 
neutron separation energies for particle-bound even-even nuclei. We used the 
SkP energy functional [16], which has a general form given by Eqs. (4) and (5). 
Our load-balancing routine, which scales the problem to 200 processors, allows 
us to perform these calculations in a single 24 wall-clock hour run on a 4 Tflop 
machine Cheetah at ORNL (1 Tflop=l x 10 12 floation-point operations/sec) 
[22]. For the details of the Skyrme-HFB deformed nuclear mass table with 
SLy4 and SkP energy density functionals, see Ref. [26]. In the following, some 
particular aspects of our Skyrme-DFT calculations are briefly discussed. 
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3. 1 Transformed Harmonic Oscillator Basis 

Going away from the beta stability valley towards particle drip lines, the Fermi 
energy becomes very small and the nucleonic densities and fields acquire large 
spatial extensions due to the coupling to the particle continuum. In this re- 
gion of weakly bound nuclei, the asymptotic behavior of nuclear densities has 
an effect on nuclear properties. Consequently, when performing calculations 
for drip-line systems, it is important to have a firm grasp on physics at large 
distances. The recently developed HFB-THO technique based on the trans- 
formed harmonic oscillator (THO) method [27,28,22] is very helpful in this 
respect: it is fast, efficient, and easy to implement. 




2 4 6 8 10 12 14 16 18 20 

r(fm) 

Fig. 3. Comparison of the neutron densities (in logarithmic scale) calculated for 
the deformed nucleus 110 Zn using coordinate-space 2D calculations (solid squares) 
with the configurational calculations based on THO (open squares) and HO (open 
circles) basis [29]. Each point corresponds to one Gauss-integration node in the 
z — p plane, and the results are plotted as functions of the distance from the origin, 
r = \J z 2 + p 2 . 

Figure 3 shows the neutron density of the deformed nucleus 110 Zn obtained in 
two configurational calculations based on expansions in the harmonic oscillator 
(HO) and THO bases [27,28,22] compared to full-fledged 2D coordinate-space 
calculations [30,29] with the box boundary conditions. Every point in the figure 
corresponds to the value of the neutron density at a given Gauss-integration 
node in the z — p plane. Since the nucleus is deformed, and there are always 
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several nodes near a sphere of the same radius r = \/ z 1 + p' 2 , there can be 
seen some scatter of points corresponding to different densities in different di- 
rections. While the significant deviation from the correct decaying behavior is 
seen in the HO results, the THO expansion agrees very well with the deformed 
coordinate-space method. Other promising techniques that can alternatively 
be used in this context are the Gaussian-expansion basis method [31] and the 
Berggren expansion method [32]. 

3.2 Regularization of the Contact Pairing Interaction 

When employing contact pairing interactions such as the density-dependent 
delta force resulting in the pairing functional (5), one has to apply a cut- 
off procedure and use a finite space of single-particle states [16]. When this 
space increases, the pairing energy diverges for any strength of the interaction; 
therefore, one has to readjust the pairing strength for each size of the single- 
particle space [17]. Such renormalization procedure is performed in the spirit of 
the effective field theory, whereupon contact interactions are used to describe 
low-energy phenomena while the coupling constants are readjusted for any 
given energy cut-off to take into account neglected high-energy effects. It has 
been shown that by carrying out renormalization for each value of the cut-off 
energy, one practically eliminates the dependence of the HFB results on the 
size of the single-particle space. 

Recently, the subject of the contact pairing force has been addressed in Refs. 
[33,34,35,36,37] suggesting the renormalization procedure can be replaced by 
a regularization scheme which removes the cut-off energy dependence of the 
pairing strength. Differences between the HFB results emerging from the pair- 
ing renormalization and pairing regularization procedures have been analyzed 
in Ref. [38] for both spherical and deformed nuclei. Figure. 4 shows differ- 
ences between the HFB-SkP results for the deformed Er nuclei obtained using 
pairing renormalization and regularization. While the regularization method 
is better theoretically motivated, it is seen that both methods give indeed very 
similar results. 



3.3 Particle Number Projection 

The advantage of the mean-field approach to the pairing problem lies in its 
simplicity that allows a straightforward interpretation in terms of pairing fields 
and deformations (pairing gaps) associated with the spontaneous breaking of 
gauge symmetry. However, in the intrinsic-system description, the particle- 
number invariance is internally broken. Therefore, to relate to experiment, 
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Fig. 4. Differences between pairing renormalization (RN) and regularization (RG) 
procedures for (a) total binding energies, (b) two neutron separation energies, and 
(d) the average neutron and proton gaps. Equilibrium quadrupole deformations are 
shown on panel (c). Calculations are performed for the chain of Er isotopes within 
the deformed HFB+THO method using SkP Skyrme parametrization and mixed 
delta pairing. 



the particle number symmetry needs to be restored. This can be done on var- 
ious levels, including the quasiparticle random phase approximation, Lipkin- 
Nogami (LN) method, the projected LN method (PLN) [39,22,21], and the 
particle-number projection before variation (PNP) [40,41,42]. 

Recently, particle-number restoration before variation has been incorporated 
for the first time into the Skyrme-DFT framework employing zero-range pair- 
ing [43]. It was demonstrated that the resulting projected HFB equations can 
be expressed in terms of local gauge-angle-dependent densities. In Ref. [43], 
results of PNP calculations have been compared with those obtained within 
LN and PLN methods. While the PLN gives results close to PNP for open- 
shell nuclei, for closed-shell nuclei it breaks down with more than one MeV 
difference in the total binding energy; see Fig. 5. This pathological behavior of 
LN and PLN methods around closed-shell nuclei can be partly cured by per- 
forming particle-number projection from neighboring open-shell systems [44]. 
This result is important in the context of large-scale microscopic mass calcu- 
lations such as those of Ref. [21]. To be on the safe side, however, it is always 
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Fig. 5. The total binding energy (with respect to a linear reference) as a function 
of the neutron number N for even-even nuclei around doubly 
100 Sn, 132 Sn calculated in LN, PLN, and PNP methods. The crosses for magic 
nuclei indicate the PLN results obtained by projecting from neighboring nuclei, as 
indicated by arrows. 

recommended to apply the complete PNP procedure around closed shells. 



4 Towards the Universal Nuclear Energy Density Functional 

Developing a nuclear density functional requires a better understanding of 
the density and gradient dependence, spin and isospin effects, and pairing, as 
well as an improved treatment of symmetry-breaking effects and many-body 
correlations. Below are summarized the areas of current theoretical activities 
in this field. 

4-1 Density and gradient dependence 

An important avenue is to enrich the density dependence of the isoscalar 
and isovector coupling constants, both in the p-h [45,46] and p-p channels 
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[47,48,49]. In particular, as the energy functional is supposed to describe those 
nuclear features that are related to collective dynamics, it seems important to 
enrich the density dependence of the effective mass in order to differentiate 
between its value in the bulk and at the Fermi surface [50,20]. 

One of the crucial challenges in microscopic theory of nuclear masses is to 
better understand salient features of the nuclear symmetry energy. The sym- 
metry energy can be extracted directly from the calculated binding energy 
of finite nuclei, after subtracting shell effects [51]. The goal is to understand 
connections between the symmetry energy and isoscalar and isovector mean 
fields, and in particular the influence of effective mass and pair correlations 
on symmetry energy versus the isospin. Such understanding will allow us to 
better determine isospin corrections to nuclear mean fields and energy density 
functionals. 

Recently, important indications on how to construct the nuclear energy func- 
tional have been obtained within the effective field theory (see, e.g., Refs. [52,53,54]). 
Even if one still has to readjust and fine-tune the parameters for a precise de- 
scription of nuclear data, one can gain important insights into the structure 
of the functional, especially the dependence of the coupling constants on nu- 
clear densities. In addition, the systematic, controlled momentum expansion 
on which the effective field theory is based offers a way to estimate theoretical 
errors (see Sec. 4.4). 



4-2 Time- odd fields 

In the self-consistent method, the average nucleonic field is obtained from 
the nucleonic density. Consequently, in a highly polarized high-spin state, the 
mean-field potential is expected to acquire appreciable time-odd components 
[55,56]. However, such terms should be present in all nuclear states with non- 
zero angular momentum, including ground states of odd-mass and odd-odd 
nuclei [57]. It is rather clear that without getting a handle on the time-odd 
fields, it will be impossible to make precise predictions for binding energies of 
most of the nuclei. 

The time-odd terms are very poorly known. An important task is to learn 
about them through an analysis of high-spin states and spin-isospin excita- 
tions. Some of the time-odd fields have been studied in Ref. [58] in the context 
of Gamow- Teller beta decays in radioactive nuclei by constraining the energy 
functional to the empirical spin-isospin Landau parameters. The coupling con- 
stants of the remaining terms can, in principle, be found by performing sys- 
tematic studies of rotating nuclei. This strategy has recently been followed in 
the Skyrme-HF analysis of high-spin terminating states [59,60]. Those fully 
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aligned states have fairly simple single-particle configurations, and they pro- 
vide an excellent testing ground for the time-odd densities and fields. 

4-3 Dynamical corrections 

The correlation term, accounting for correlations going beyond the simple 
product state, is an integral part of the DFT. Since nuclei are self-bound sys- 
tems, many-body correlations due to spontaneous symmetry-breaking effects 
are of particular importance. A large part of those correlations can indeed be 
included by considering symmetry-breaking product states. Within the mean- 
field approach, one can understand many physical observables by directly em- 
ploying broken-symmetry states; however, for finite systems, a quantitative 
description often does require symmetry restoration. For this purpose, one 
can apply a variety of theoretical techniques, in particular projection methods, 
the generator coordinate method (GCM), the random phase approximation 
(RPA), and various approximations performed on top of self-consistent mean 
fields [61,62,63]. 

In this context, it is important to recall that the realistic energy density func- 
tional does not have to be related to any given effective Hamiltonian. This 
creates a problem if a symmetry is spontaneously broken. While the projec- 
tion can be carried out in a straightforward manner for energy functionals 
that are related to a Hamiltonian, the restoration of spontaneously broken 
symmetries of a general density functional still poses a conceptional dilemma 
that needs to be properly addressed [64,65,66]. 

Since the correlation term is a part of the functional, it should be treated as 
such during the variational procedure and during the fitting process in which 
the functional's coupling constants are determined. So far, perhaps with the 
exception of the center-of-mass term (see Sec. 4.3.1 below), such an ambitious 
program has not been carried out. In the near future, one hopes to work out ap- 
proximate expressions for the correlation term that would capture the essence 
of results of microscopic calculations performed on top of self-consistent mean 
fields. In this way, the hope is to develop the tractable parametrization of 
the correlation energy in terms of local densities that would allow an explicit 
inclusion of dynamical effects into the energy functional. 

4-3.1 Center-of-mass correction 

The center-of-mass (cm.) correction, due to the violation of the translational 
invariance, is always included in calculations, but its practical implementations 
differ from functional to functional [57,8]. For some functionals the treatment 
is fully variational; for some others the cm. term is computed following the 
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HFB procedure; for some functionals a simple one-body approximation is used. 
These apparently technical differences do matter as the actual form of the cm. 
correction has a significant impact on the surface properties [57]. 

A good example nicely illustrating the above point has recently been discussed 
in Ref. [51]: for the two functionals, SLy4 and SLy6, which were fitted with 
precisely the same strategy but differ in their treatment of cm. correction, the 
surface energy coefficient differs by as much as 0.7 MeV. While the two-body 
(albeit perturbative) treatment of the cm. correction does not reduce the 
overall rms error of the fit to nuclear masses [20], it certainly has a significant 
impact on binding energies of highly deformed configurations (such as fission 
isomers), fission barriers, and fission trajectories. 

4-3.2 Particle-number and isospin corrections 

As already discussed in Sec. 3.3, efficient numerical codes that allow for large- 
scale, self-consistent variational calculations after projecting onto a good par- 
ticle number have been developed [43]. The particle-number conserving HFB 
equations [40,41] with Skyrme functionals can be simply obtained from the 
standard Skyrme-HFB equations in coordinate space by replacing the intrin- 
sic densities and currents by their gauge-angle dependent counterparts. Using 
the VAP method, one can properly describe transitions between normal and 
superconducting phases in finite systems, which are inherent in (semi) magic 
nuclei. 

As mentioned above, the restoration of broken symmetries in the framework of 
DFT causes a number of questions, mainly related to the density dependence 
of the underlying interaction and to different treatment of particle-hole and 
particle-particle channels [42,67]. These questions are a matter of ongoing 
intensive research [65,66]. 

Related to the particle-number symmetry, but different in origin and treat- 
ment, is the question of the spontaneous isospin breaking. The isospin-breaking 
correction is of particular importance around the iV~Z line. The isoscalar 
pairing is believed to contribute to the additional binding of N=Z nuclei, the 
so-called Wigner energy [68]. However, basic questions regarding the collectiv- 
ity of such a phase still remain unanswered, and should be part of the future 
scientific agenda. 

Apart from the presence of charge-dependent terms in the functional, such as 
the Coulomb term, the isospin symmetry is broken by the quasiparticle mean 
field (the generalized product wave function is not an eigenstate of isospin). 
Several techniques have been developed to restore isospin (see the discussion 
in Refs. [10,69] and references quoted therein). It is fair to say, however, that 
in spite of many attempts to extend the quasiparticle approach to incorporate 
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the effect of proton-neutron correlations, no symmetry-unrestricted mean-field 
calculations of proton-neutron pairing, based on realistic effective interaction 
and the isospin-conserving formalism, have been carried out so far. 

4-3.3 Rotational and vibrational zero-energy corrections 

The rotational-vibrational correlations are important aspects of nuclear col- 
lective dynamics; they also contribute to nuclear binding through quantum 
zero-point corrections. To estimate the magnitude of the rotational-vibrational 
corrections, one usually applies RPA [70], GCM [63], or the Gaussian overlap 
approximation to GCM [71,72,73,74]. 

Regardless of the approach used, a key point is the choice of collective sub- 
space. In the case of GCM and related methods, the collective manifold is 
determined by the set of external fields associated with the collective motion 
of the system. In most practical applications, one considers five quadrupole 
degrees of freedom that give rise to nuclear rotations and quadrupole vibra- 
tions, octupole deformations, and pairing vibrations [75,70]. An important 
step towards the microscopic description of correlation energies are the recent 
large-scale benchmark calculations of ground-state quadrupole correlations of 
binding energies for all even-even nuclei, from 16 up to the superheavy sys- 
tems [63]. 

4-4 Fitting Strategy and Error Analysis 

One of the still-unsolved questions is an appropriate selection of experimental 
data that would allow for a more-or-less unique determination of the coupling 
constants defining the energy functional. To this end, one usually uses certain 
constraints obtained by extrapolating nuclear data to an infinite system and 
selected data for finite nuclei. The sensitivity of the final fit to the choice of 
this data set leads to a plethora of parameterizations currently available in 
the literature. 

Most of the currently used density functionals correctly reproduce generic 
trends in nuclear masses - as selected masses are usually considered in the 
data set - but their descriptions of other quantities vary. Moreover, they of- 
ten significantly differ in parameters or coupling constants [8]. This suggests 
that yet-unresolved correlations may exist between these parameters, and only 
certain combinations thereof are important [76,51]. Such correlations would 
explain the fact that widely different parameterizations lead to fairly similar 
results. 

The present stage of theory requires constructing new energy density func- 
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tionals supplemented by a complete error and covariance analysis. It is not 
sufficient to "predict" properties of exotic nuclei by extrapolating properties of 
those measured in experiment. One must also quantitatively determine errors 
related to such an extrapolation. Moreover, for experimental work it is essen- 
tial that an improvement gained by measuring one or two more isotopes be 
quantitatively known. From a theoretical perspective, one must also know the 
confidence level with which the parameters of the functional are determined. 
An analysis of this type constitutes a standard approach in other domains 
of physics, but they are seldom performed in theoretical nuclear structure 
research. 



5 Summary 

This paper discusses the status, advances, open problems, and perspectives 
in the area of large-scale microscopic nuclear mass calculations. This field of 
research is past phenomenological approaches that gave us a very good under- 
standing of general features and trends, but lacked fundamental derivations 
and had limited predictive power. At present, the focus is on microscopic de- 
scriptions of nuclei whereupon they are treated as finite quantum objects built 
of (quasi) nucleons. Nuclear ground states and masses are in this approach de- 
termined by basic fields, which are the particle and spin densities along with 
their derivatives and gradients up to the second order in relative momenta. 
These fields interact in such a way that the total energy of a given system 
is a functional of densities, defined and understood in the general framework 
of the Kohn-Sham theory. The determination of such a universal functional, 
along with all the dynamic corrections required by data, is the main purpose 
of current investigations. In this endeavor, we strive not only to have a reliable 
theoretical tool to calculate properties of very exotic systems that will not be 
soon accessible in experiment, but also wish to have a spectroscopic quality 
description of well-known systems, in which very precise data do exist now, 
and can be used as a rich source for determination of theoretical parameters. 
Such a program of research should npt only be rooted in the fundamentals 
of low-energy QCD methods and ideas, but also, by definition, must rely on 
experiment for elements that cannot be derived from first principles. It is a 
vast and ambitious program presently under way. 
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